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Abstract — An exact expression for the symbol error probability 
of the four-dimensional 24-cell in Gaussian noise is derived. 
Corresponding expressions for other regular convex polytopes 
are summarized. Numerically stable versions of these error 
probabilities are also obtained. 



I. Introduction 

THE 24-CELL is a regular convex polytope (RCP) with 
24 facets and 24 vertices in four dimensions. Its vertices 
were considered as a signal constellation for digital modulation 
already in 1977 Q |. The symbol error probability (SEP) of this 
modulation format was approximated by a union bound in [T) 
and the bit error probability was estimated by simulation in Q. 
The constellation was considered for satellite communications 
in El, [3] and has recently received renewed interest in fiber- 
optical communications 0), @. This new interest was moti- 
vated by the rapid progress in coherent detection, which makes 
it possible to transmit data in both polarizations simultaneously 
and thus to employ four-dimensional constellations (6). In 
this paper, we calculate its exact SEP over the additive white 
Gaussian noise (AWGN) channel, which has not been done 
before. 

We put the result into context by briefly summarizing what 
is known about signal constellations constructed from the 
vertices of other RCPs in n dimensions. An RCP is a convex 
polytope with the maximum degree of symmetry; its facets and 
vertices all look the same, and the facets are in turn RCPs. 
Therefore, all symbols have the same energy and the same 
conditional error probability, which is tractable at least from 
a theoretical viewpoint Q. A drawback is that the number 
of vertices M of an RCP is usually not a power of two, 
but such constellations can nevertheless be used in digital 
communications if a block of bits is mapped to more than 
one symbol JS), (9), in coded modulation |10|, or to reduce 
the peak power in orthogonal frequency-division modulation 
systems ifTTI . Furthermore, in a data communication scenario, 
some levels beyond a power of two are useful for framing 
and control purposes; indeed, 3- and 5-level modulations are 
already applied in standards for Fast Ethernet and Gigabit 
Ethernet fT2] pp. 285-289]. 

The SEP of RCPs can be approximated by sphere packing 
0, random coding J7), simplex packing JT3J , and union 
bounds lfT4l Sec. 3.2.7]. To obtain an exact expression for 
the SEP, one needs to integrate an n-dimensional Gaussian 
density, centered at an RCP vertex, over its Voronoi cell. For 
a given RCP, the Voronoi cells are unbounded n-dimensional 
pyramids with apex at the origin, whose (n — 1) -dimensional 
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base cells are the facets of the dual polytope. Thus, the exact 
SEP is in general an n-dimensional integral, but as we shall 
see, it can in most cases be simplfied into a one-dimensional 
integral, thanks to the inherent symmetries of RCPs. 

There exist infinitely many RCPs in dimension n = 2, 
five for n = 3, six for n = 4, and three for n > 5 lfT31 
p. 412]. For most of these constellations, the corresponding 
SEP is known, but at high signal-to-noise ratios, the published 
expressions often consist of the difference between two almost 
equal entities, which makes them numerically unstable. We 
therefore derive equivalent expressions that lend themselves 
better to numerical evaluation. 

II. Power and spectral efficiency 

We define the signal-to-noise ratios 7 = E s /Nq and 7b = 
E\,/No, where E s is the symbol energy, E\, = E s /\og 2 M is 
the bit energy, M is the number of levels, and Nq/2 is the 
double-sided spectral density of the AWGN. For modulation 
using the vertices of any RCP, the SEP P can be upperbounded 
by the union bound 



(1) 



where Q(x) = J°° f(x)dx, f(x) is the normalized Gaussian 
probability density function f(x) = (l/\/27r) cxp(— x 2 /2), d 
is the minimum Euclidean distance between pairs of vertices, 
and A is the number of such pairs involving a given vertexQ 
The bound is an accurate approximation at medium to high 7 
and it is asymptotically exact. We define the power efficiency 
as Gb = d 2 /(4i?b), which is the reduction of bit energy E\, 
that a particular constellation admits, compared with binary 
phase-shift keying (BPSK) modulation, to attain the same, 
asymptotically low, SEP. Furthermore, the spectral efficiency 
S = log 2 M/n is the reduction of bandwidth that a particular 
constellation admits, compared with BPSK, with the same bit 
rate and pulse shape. 

The parameters M, A, and d of all RCPs are given in Tab. U 
adapted from [Tgl pp. 292-2950, along with G b and S. 

III. Generic RCPs 

Three types of generic RCPs exist for all n > 2. First, 
the n-cube can be defined as the vertex set (±1,...,±1), 
with all 2™ sign combinations. When used for modulation, 
this constellation results in independent BPSK modulation in 
all n dimensions. It specializes into a square and cube for 

'The inequality does however not hold for irregular polytopes, for which 
in general not only the closest neighbors contribute to the Voronoi cell. 

2 The parameters were calculated from jVrj, N\, and <j> in (IfSjj as M = No, 
A = 2Ni/N , and d/^/El = 2 sin <f>. 
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n = 2 and n = 3, resp. Second, the n-crosspolytope (or 
orthoplex), whose 2n vertices can be taken as all permutations 
of (±1, 0, . . . , 0), is the geometric dual of the n-cube and 
corresponds to biorthogonal modulation. It is also a square 
for n — 2 but an octahedron for n = 3. The third type is 
the self-dual n-simplex, which is an equilateral triangle for 
n = 2 and a regular tetrahedron for n = 3. Its n + 1 vertices 
can be taken as (1 — a, . . . , 1 — a) and all permutations of 
(a — n, L . . . , 1), where a = y/n + 1 and all vectors have 
length n0 

The SEP of these three RCPs are given by the following 
theorem, where two expressions are given for each SEP; 
one that is suitable for analysis and another for numerical 
evaluation. Expressions related to our © and but less 
reliable at very high 7, were given in ifTTl . 

Theorem 1: The SEPs of the n-cube P c , n-crosspolytope 
P cp , and n-simplex P s are 



TABLE I 
Parameters of All RCPs. 




\{-Q{x)fdx 



i=0 



Proof: Transmission over the n-cube is equivalent to 
transmitting a block of n bits over independent BPSK chan- 
nels, which yields ©. Expressions © and © were derived 
from (Til eqs. (4.102) and (4.116)] via the substitutions 
q = x/y/2 — yfi and q = x/\/2 — + 1/n), resp. 

Expressions ||3}, (0, and (0 follow after expanding (1— Q{-)) n 
by the binomial theorem and integrating out the constant term. 

□ 

Other geometric properties of the three generic RCPs were 
calculated in ifTHl. 



3 To our best knowledge, this is the first time that an explicit n-dimensional 
coordinate representation is presented for the n-simplex. (a = — yjn + 1 
gives another, equally valid, representation.) The standard representation, 
which is in n + 1 dimensions, is less useful for transmission. 
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IV. Two-dimensional RCPs 

In two dimensions, RCPs are polygons, which exist with 
an arbitrary number of vertices M > 3. Using such a 
constellation for transmission corresponds to Af-ary phase 
shift keying (Af-PSK). The following expression, which was 
given in lfl9ll . is both simple and stable. 

Theorem 2: The SEP of the M -polygon Ppsk is 



Pp: 
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exp 
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The SEP for selected values of M is plotted in Fig. Q] (a). 
As expected, M — 4 has the same asymptotic performance 
as BPSK. The only Af-PSK format with a power efficiency 
higher than dB is 3-PSK, as shown already in [20|. 

V. Three-dimensional RCPs 

In three dimensions, there are two RCPs in addition to 
the above-mentioned cube, octahedron, and tetrahedron. These 
are the icosahedron and the dodecahedron, with 12 and 20 
vertices, respectively, which are duals of each other. Their 
exact SEPs are not known, but their union bounds follow 
from ([TJ and Table U The results are shown in Fig. Q] (b). 
The tetrahedron (3-simplex) has the highest power efficiency 
of the three-dimensional RCPs Q. 

VI. Four-dimensional RCPs 

Coxeter remarked that a "peculiarity of four-dimensional 
space is the occurrence of the 24-cell. . . , having no analogue 
above or below" lfl6l p. 289]. The 24-cell is a four-dimensional 
self-dual RCP with 24 vertices. It is the only RCP, in any 
dimension, that is both more power efficient and more spec- 
trally efficient than the n-cube (BPSK modulation). Curiously 
enough, the power and spectral efficiencies of the 24-cell are 
equal, d 2 /(AE b ) = log 2 M/n = (3 + log 2 3)/4 = 1.15. 

The 24-cell can be described in several ways depending 
on its orientation. The vertices can be taken as the union of 
(±1,±1,±1,±1), with all possible combinations of signs, and 
11(2, 0, 0, 0), where II(-) means the set of all permutations of 
the given elements. In other words, the 24-cell is the convex 
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hull of the union of the 4-cube and the 4-crosspolytope, scaled 
to the same circumscribed radius and properly aligned. This 
description, although geometrically appealing, does not seem 
to easily admit an exact computation of the SEP. 

Rotating the constellation, the vertices of the 24-cell can be 
described as the set 

x = n(±i, ±i,o,o), 

with all possible permutations and signs. With this coordinate 
representation, we are able to calculate the SEP of the 24-cell. 
The result is again presented in the form of two expressions, 
where the first is analytically simpler and the second numeri- 
cally useful even at high 7 values. 
Theorem 3: The SEP of the 24-cell is 



With these definitions, 



1 



exp 



{l-2Q(x)) 2 Q(x-^)dx 



Q(V7)(2-Q(V7)) 



exp 



(8) 

{x - V7) 2 



Q(x)(l-Q(x))Q(x~^)dx. 



(9) 



Proof: A signal vector X\/E s /2, where X is ran- 
domly taken from X with equal probabilities, is transmitted 
over a discrete-time AWGN channel with variance No/ 2. 
The received vector is X^jEj2 + Zy/N Q /2, where Z = 
(Zi, Z 2 , Z 3 , Z4) is a vector of independent, zero-mean, unit- 
variance, Gaussian random variables. To simplify the analysis, 
the received vector is rescaled by a factor ^jE s /2. Thus, the 
maximum likelihood detector shall find the vector in X being 
closest to X + Z/i/j in the Euclidean sense. 

Because the signal vectors are equally probable and sym- 
metrically located, we assume, without loss of generality, that 
the transmitted vector is x 4 = (1,1,0,0). The SEP can now 
be expressed as 



P 24 = Pr{xi + — Z i n*(xi)}, 
where D,x{u) for any it G R 4 is the Voronoi celfl 



(10) 



Q x {u) = {y G 



\y ■ 



<\\y-x\\,xeX}. (11) 



In order to integrate the four-dimensional noise over the 
Voronoi cell, we need to find a compact coordinate represen- 
tation for this cell. To this end, we first partition X into the 
subsets 

X x 4 {(1,1,0,0)}, 

Af 2 4{(0,l,0,±l),(0,l,±l,0),(l,0,0,±l),(l,0,±l,0)}, 
X 3 = {(-1,1, 0,0), (1,-1, 0,0)}, 

,*4 = {(0,0,±1,±1)}, 

# 5 ^{(0,-l,0,±l), (0,-1, ±1,0), 

(-1,0,0,±1),(-1,0,±1,0)}, 
X 6 4 {(-1,-1,0,0)}. 

4 In our notation, the definition is valid also when u £ X, which will be 
utilized in (U\-(Ti\. 



X = \Jx. 



i=l 
6 



Flx{u) = (^1 ftxiiu), for any u E 



(12) 



Applying (fTTT i to the subset X 2 yields 

Qx 2 {x-l) 

= {yeR 4 I Wy-x.W 2 < \\y - x\\ 2 ,x e X 2 } 

= {y e R 4 1 (y, x 1 -x)>0,xe X 2 } 

= {y e k 4 I yi > M,2/i > Is/3 1, 2/2 > Is/4], s/a > I2/3] > 

= {y G R 4 I mm{ yi ,y 2 } > max{| 2/3 |, \y 4 \}}, (13) 

where y = (1/1,2/2,2/3,2/4) and (•, •) denotes the inner product. 
The Voronoi cells of X\ with respect to the subsets X3, . . . , Xq 
can be calculated by the same method. Omitting the details, 
the results are 



iW*i) - {y E M 4 I mm{ yi ,y 2 } > 0} 



(14) 



iW*i) = {y e R 4 1 2/1 + 2/2 > I2/3I + I2/4I} (15) 

toxaixi) = {y e M 4 I 2/1 +2/2 +min{2/i,2/ 2 } 

>max{|2/ 3 |,|2/ 4 |}} (16) 
O* 6 (zi)={yeR 4 |2/i+2/2>0}. (17) 

Comparing (\3[ with (fl4l yields Slx 2 ( x i) C Q^ 3 (xi), be- 
cause min{ 2/1,2/2} > max{|2/ 3 |, I2/4I} => min{y 1 ,y 2 } > 0. 
Similarly, comparing (fT3l l with each of (fT5ll-(fT7]i shows that 
£lx 2 ( x i) Q ^Xi( x i) f° r i — 4,5,6. Now using these results 
in (fT2l yields 

Qx(xi) = tl X2 (xi). (18) 

We can now finish the SEP calculation. From ( TTOb , dT3l . 
and CD, 

Pu = 1 - Pr jxi + -^=Z € fi* 2 (xi) 

= l-Pr(niin(l + ^.l + ^l>niax(^,^l 

= l-Pr{V7 + min{Zi,Z 2 } > max{|Z 3 |, \Z 4 \}} 
= l-Pr{Zi < Z 2 , V7 + min{Zi,Z 2 } > max{|Z 3 |, |Z 4 |}} 
-Pr{Zi > Z 3 ,y/i + mm{Zi,Z 2 } > max{|Z 3 |, |Z 4 |}}. 
By symmetry, the last two probabilities are equal; thus, 

P 24 = 1 - 2Pr{Zj < Z 2 , V7 + min{Zi, Z 2 } > max{|Z 3 |, |Z 4 |}} 
= l-2Pr{Zi < Z 2 ,V^+Z! > \Z 3 \,y/7 + Z x > \Z 4 \}. 

Finally, we marginalize the probability over Z\ = z\ to obtain 



Pu = 1 - 2 



•Pr{Z 2 > zi,|Z 3 | < yfi + zx,\Z A \ < y/j + zijdz! 

/oc 
/(«i)Q(*i)(l - 2Q(^7 + Zi)) 2 ^! (19) 

/» OC 

= 1 - 2 / /(z - V7)Q(^ - V7)(l ~ 2Q(z)) 2 dz, 
Jo 
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Fig. 1. The SEP of RCPs in various dimensions n. Except for the icosahedron, dodecahedron, 600-cell, and 120-cell, whose SEPs are approximated by the 
union bound, all SEPs are exact. The SEP of BPSK is shown as a reference in all diagrams (dashed line). 



which completes the proof of (H). 

To prove (O, we first note that for any x E R, 

Q 2 {x) = Pt{Z 1 >x,Z 2 > x} 
= 2Pt{Z 2 >Zi>x) 

= 2 / f{z 1 )Q{z 1 )dz 1 . (20) 

J X 

Expanding < [T9l > with the help of (l20t yields 

/oo 
f(zi)Q(zi) 

• (1 - 4Q(V7 + zi) + 4g 2 (V7 + zi))d Zl 

/OO 
f{zx)Q{z 1 ) 

• (4Q(V7 + *i) - 4Q 2 (V7 + zx))dzi 

/CO 
f(zMzt) 
-Vv 

■ QiVl + z^il - Q(V7 + zi))dz! 
= Q(Vt)(2-Q(V7)) 

+ 8 / /(z - v^)Q(z - Vt)Q(2)(1 - Q(z))dz, 
Jo 

which completes the proof of (O. □ 
The SEP of all four-dimensional RCPs is shown in 
Fig. [T] (c). It can be observed that n = 4 is the lowest 
dimension for which the crosspolytope has the highest power 
efficiency of all RCPs. In addition to the 24-cell and the three 
generic poly topes, there are RCPs with 120 and 600 vertices in 



four dimensions. Miyazaki made an artistic effort to visualize 
these polytopes ||2T1 Ch. 14], which are each other's duals. 
Their power efficiencies are low, similarly to the two largest 
RCPs in three dimensions. 



VII. Higher-dimensional RCPs 

In dimensions n > 5, there are no RCPs apart from the three 
generic types. For medium and high 7b, the crosspolytope is 
always better than the simplex, and the simplex is better than 
the cube. Their SEP performances are shown in Fig. G](d) for 
the five-dimensional case. 

As the dimension increases, the power efficiency gap be- 
tween the simplex and the cube increases monotonically. 
However, the gap between the crosspolytope and the simplex, 
which based on Tab. U and the definition of Gb can be 
expressed as 

Gb,c P = nlog 2 (2n) 
G b ,s (n + l)log 2 (n + l)' 

increases until it reaches a maximum for n = 24, where 
the crosspolytope is 0.62 dB better. Thereafter the gap de- 
creases again and approaches dB in very high dimensions. 
Extending n to real numbers, the maximum of ( 1211 occurs 
for n 24.066, which is quite close to 24. Whether this is 
a pure coincidence or more deeply related to other unique 
properties of 24-dimensional geometry, such as the existence 
of the Leech lattice, is not known. 
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